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ABSTRACT 



We present results of our survey for planetary transits in the field of NGC 6940. 
We think nearly all of our observed stars are field stars. We have obtained high pre- 
' cision (~3-10 millimags at the bright end) photometric observations of ^50,000 stars 

spanning 18 nights in an attempt to identify low amplitude and short period transit 
q ■ events. We have used a matched filter analysis to identify 14 stars that show multi- 

$H ' pie events, and four stars that show single transits. Of these 18 candidates, we have 

, identified two that should be further researched. However, none of the candidates are 

convincing hot Jupiters. 

Key words: methods: data analysis - stars: variables - open clusters and associations: 
individual (NGC 6940) - planetary systems 
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1 INTRODUCTION 

ICharbonneau et alJ l|200(t) opened a new chapter in the sci- 
ence of extrasolar planets when they recorded the first tran- 
sit of a planet around its parent star. The transit produced 
a 1.5% dip in the star's light. Until then, the only evidence 
of planets around main sequence stars had been radial ve- 
locity measurements of stellar reflex motions. Though the 
RV method has been the most successful method of find- 
ing planets heretofore, the transit method of searching for 
planets is complementary, because it provides different infor- 
mation than RV. Measuring a transiting planet can provide 
the actual mass of the planet by determining the orbital in- 
clination of the system and provide the radius of the planet. 
Also, in some situations, transiting planets c an be probed 
for a tmospheric spectra, as with HD 209458 jBrown et alJ 
l200ll) . Finally, the transit method can find planets to kilo- 
parsec distances, much farther than RV. 

However, the strength of the transiting method of dis- 
covery, that it shows us the orbital inclination, is also 
its weakness, because that orbital inclination must be 
close to 90 degrees for us to see the transit. Radial ve- 



locity measurements have shown that approximately 1- 
2% of Sun-like stars in the solar neighborhood have hot 
Jupi ters, giant planets with orb ital distances of 0.035-0.4 
AU |L mcwcavcr fc GretheJ2003Th Assuming that orbital in- 
clinations are random, approximately 10% of stars with hot 
Jupiters should have transits visible to us. Therefore, ap- 
proximately one of a thousand Sun-like stars should show 
an eclipse, if the stars we observe have the same planetary 
abundance as the solar neighborhood. 

Ijanesl il996l) suggested that open clusters would be 
good fields in which to look for planetary transits. Open 
clusters contain hundreds of stars of a similar distance and 
metallicity. The field is crowded enough to be able to observe 
a sufficient number of stars, but not crowded enough to make 
reduction exceptionally difficult. The high number of stars is 
essential, since perhaps only one in a thousand stars will ex- 
hibit the characteristic dip (a shallow flat-bottomed eclipse) 
of a planet transiting the parent star. Unfortunately, though 
this is the reason we observed in the direction of NGC 6940, 
we don't think we have observed any significant number of 
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Table 1. Parameters of open cluster NGC 6940 

RA (J2000.0) 20 34 26 

Dec(J2000.0) +28 17 00 

I 69.90 

b -7.17 

Distance(pc) 770 

Distance modulus (mag) 10.10 

Age(Zogio) 8.858 

Age(Gyr) 0.72 

[Fe/H] +0.01 

E(B-V) 0.214 



cluster members. We describe the reasons for this in more 
detail in section 2.5 below. 

We present results from a deep search for planetary 
transits in the field of NGC 6940. We describe the obser- 
vation and data reduction methods used in order to extract 
light curves for each of these stars. We show that using these 
methods we can achieve the accuracy necessary to detect 
planetary transits of a Jupiter-radius object. We describe 
our transit finding algorithm and show with simulations that 
we can recover injected transits using that algorithm. Fi- 
nally, we describe several transit candidates: 14 stars that 
show multiple low amplitude short duration events and four 
stars that show single events. We have rejected all but two 
as poor transit candidates, and recommend them for further 
study. 



2 OBSERVATIONS AND REDUCTION 
2.1 Observation 

Observations were taken over June and July of 1999 using 
the 2.5 metre Isaac Newton Telescope at La Palma, Canary 
Islands. Usable observations were taken on 18 nights be- 
tween 22-30 June and 22-31 July. Images were taken with 
the Wide Field Camera, a mosaic consisting of four 2048 x 
4096 pixel EEV CCDs, mounted at the prime focus of the 
INT. The mosaic created a 0.29 square degree field of view 
with 0.33 arcsec per pixel (see Fig. 0. 

Three open clusters were o bserved in rotation during 
the observing run, N GC 6819 JStreet et al.1 l200St . NGC 
6940, and NGC 7789 feramich et alj|2004h . This paper re- 
ports on the analysis of NGC 6940 (see Table 1) observa- 
tions. Each image was exposed for 300 seconds, taken in 
pairs to help remove/identify cosmic rays. This resulted in 
approximately 2 observations per hour per cluster. We ob- 
tained 251, 278, 267 and 249 usable frames of NGC 6940 for 
each of the four CCDs, respectively. The observing routine 
was designed to maximize the number of stars observed, in 
order to maximize the possibility of a transit detection. The 
300 second exposure setting was mainly in order to capture 
enough cluster member stars of NGC 6819 and NGC 7789, 
which are 1900 and 2400 parsecs distant, respectively. This 
setting has caused some minor problems with the observa- 
tion of NGC 6940, discussed below in the section on colours. 
In retrospect, a shorter exposure time would have been bet- 
ter for NGC 6940, to avoid saturating cluster stars at 770 
pc. 




Figure 1. CCD Mosaic of NGC 6940. 



2.2 Data Reduction 

After standard CCD processing, the individual science 
frames were reduced with differential im age analysis, based 
on code developed bv lBond et alJ l 200 if). The process i s de- 



scribed in more detail by Bramich ( Bramich et all2004l) and 
summarized here. 

We used an automated script and IRAF tools to build a 
3-sigma clipped mean masterbias and 3-sigma clipped mean 
masterflat frame. From each of the science frames we then 
subtracted the masterbias frame and divided the master- 
flat. For the reduction procedure, we considered each of the 
CCDs separately. However, unlike Bramich, we considered 
all the observations as one run, over June and July 1999, 
instead of considering them as separate runs. 

Following the standard processing, we reduced the pho- 
tometry on the science frames using differentia l image anal- 
ysis (DIA) llAlard fc Luptorjll998tlAlar dll2000D . Our imple- 
mentation o f DIA code was written for the MOA project 
(Bon d~t alJ l200lf) . All of the processes are automated 
into scripts which call on C code developed by Bond and 
Bramich. 

Differential Image Analysis (DIA) is excellent for accu- 
rately measuring variable stars within a somewhat crowded 
field. The idea of differential image analysis is that constant 
stars are removed from the observations, leaving only those 
stars in which we are interested, because they contain vari- 
ability induced possibly by a transiting planet. We first used 
a script to build a reference frame tha t is a combin ation of 
the best seeing frames in the entire run. lAlardl (l2000l) showed 
that using several good seeing science frames generated bet- 
ter results than just using one, best frame as the reference. 

We subtracted this reference frame from each of the 
science frames to create residual images. In order for the 
subtraction to be successful, we had to convolve the refer- 
ence frame to the same seeing as each of the science frames. 
The science frames I(x,y) are related to the reference frame 
R(x, y) with the convolution equation: 
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I(x, y) = K(u, v,x,y)® R(x, y) + B(x, y) (1) 



where K(u,v,x,y) is the convolution kernel and B(x,y) 
represents the sky background. Thus, the residual images 
should have only random noise at the positions of constant 
stars, while the variable stars will create a dark or light 
spot on the residual, depending on whether or not the star 
was dimmer or brighter (relative to the reference frame) in 
the working image. This method generally performs much 
better than PSF fitti ng, particularly with blended stars 
flAlard fc LuD tonl ll99£h . 

Finally, we measure the flux on the residual images us- 
ing an optimal PSF scaling at the position of each star. 
Stars have already been identified using PSF fitting on the 
reference image, using IRAF's DAOPhot package. 



2.3 Photometric Precision 

We find that with the above processing, we can achieve an 
rms scatter of 0.004 - 0.006 mag at the bright end of our 
observations; good enough to detect planetary transits (see 
Fig-HJ- However, only a very small number of our stars have 
precision near this limit. Only ~4400 stars of the ~50,000 
have rms scatter better than 1%. 

Our instrumental magnitude saturation limit for each 
of our CCDs was approximately 17. Beyond that limit, sat- 
urated stars, bad columns, and CCD defects were identified 
as stars. We also see that CCD three (Fig. 2 c) has a much 
tighter curve than the other three CCDs. This is because we 
were able to combine 12 best seeing frames in order to make 
the reference frame for CCD three. The constituent frames 
of the reference frame need to be roughly sequential, or at 
least occur on the same night, and only CCD three had such 
a run of sequential, good seeing frames, without defects. The 
other CCDs only had four to six sequential frames with good 
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(a) CCD 1 (b) CCD 2 




(a) CCD 3 (b) CCD 4 

Figure 3. Colour magnitude diagrams. Colour magnitude diagrams for each of the four CCDs, with the colours and magnitudes converted 
to standard values. The highlighted stars are our transit candidates. The line represents a theoretical main sequence for a cluster 770 
parsecs away, but only K and M stars would be represented by the line. The main dark ridge in each of the graphs, with R — I colour 
indices of 0.5-0.75 are spectral type K3-K8. 



seeing (most had output errors) . This created slightly worse 
reference frames on the other three CCDs and stars with 
more scatter in the measure of precision. 



2.4 Colour Data 

We found ~350 photometric standard stars for the field 
of NGC 6940 from the Canadian Astronomy Data Centre 
JStetsonll2000h . Of these standard stars, we were able to 
use ~240 stars to calibrate the observations (~110 of the 
standard stars were saturated in our data). To change our 
observations from the instrumental CCD magnitudes into 
standard R and / magnitudes, we made a linear regression 
to put CCD one into the standard observations, then cor- 
rected each of the other CCDs to conform roughly to CCD 
one's values. 

We began by computing a linear regression between 



the instrumental r and i values and the standard (Johnson- 
Cousins) observed R and I values of our 240 stars: 

Rjc = rccD x 0.977 + 0.192 (2) 

and 

lie = «ccd x 0.985 - 0.702 (3) 

Unfortunately, the photometric standard stars observed 
were in the center of the cluster, so they only appear on 
CCD one. Thus, offsets were inferred for the remaining three 
CCDs by assuming that the mean magnitude in r and the 
colour r — i of all the stars (to magnitude 20, when we have 
large errors) would be approximately equal. We found that 
the following offsets correct the biases of the other CCDs: 

(r - i)ccD2 = (R-I) + 0.064 (4) 

r C CD2 =R- 0.105 (5) 
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Figure 4. Photometric precision vs stellar radii. The lines show 
the transit amplitude that would occur with planets 0.5, 1.0, and 
2.0 Rj up . 



(r - i) CC D3 = (R -I)- 0.338 (6) 

r C CD3 =R- 0.177 (7) 

(r - i)ccD4 = (R -I)- 0.073 (8) 

r C CD4 = R- 0.148 (9) 

2.5 Colour Magnitude 

Using the calibrations above, we computed the R— I colour 
index for each star and produced a colour-magnitude dia- 
gram (Figure 3). We were unable to find a significant main 
sequence in the observations of NGC 6940. The 300 s ex- 
posures have saturated the members of our cluster, which 
is only ~770 pc distant, as opposed to the much larger dis- 
tances of the other clusters (~1900 and ~2400 pc). Only the 
K and M cluster members of NGC 6940 were faint enough 
to be observed, but we believe that not enough of these 
have been observed to consider our stars part of the cluster. 
We assume that all our data only refers to field stars. The 
Besangon model for our direction of the galaxy and our ob- 
servation limits in R estimates that 94.5 % of our observabl e 
stars should be K0 spectral type or later (iRobin et all2003l) . 
This correlates very well with the observed R — I colour in- 
dexes for our stars, which suggest that 94.3% of our stars 
are of K0 spectral type or later. 

2.6 Stellar Radii 

We used the calibrated R — I index to estimate the stellar 
radius for each of the stars in our data set. We did this by 
interpolating be tween stan dard values of stellar radius and 
standard R — I JCoxll2000t to arrive at this polynomial: 

R,/R Q = 1.333-1.548(i?-7)+1.131(i?-/) 2 -0.3501(7?-/) 3 (10) 

Judging from the calibrated R — I index, nearly all our 
main sequence stars are K and M type. We can determine 



roughly the stars that have sufficient precision by compar- 
ing the rms of the star with the depth of a theoretical tran- 
sit of a Jupiter-sized object. Figure 4 shows the scatter of 
each of our stars compared to the stellar radii. The lines 
represent 0.5 Rj up , 1 Rjup, and 2Rj up transits in front of 
stars with the appropriate stellar radii. Almost none of our 
stars have precision good enough to view the transit of a 0.5 
Rjup planet, however, about 19% of our stars have enough 
precision to measure a one Rj up transit, while nearly 56% 
have precision to measure a two Rj up transit. We have not 
provided a rigorous treatment of extinction and reddening, 
which will affect the computed size of the stars, because we 
are observing field stars, with varying distances. 



3 TRANSIT DETECTION ALGORITHM 

The final step in our data reduction is the search for plan- 
etary transits from the stellar light curves. We used a 
matched filter algorithm which compares theoretical transit 
lightcurves with the observed lightcurves from our ~50,000 
reduced stars. 

This search uses a truncated cosine approximation with 
four parameters: period, duration, depth, and the time of 
transit midpoint. We first used a period sweep from 1.5 d to 
7 d with a fixed transit duration of 3 h. The stars with mul- 
tiple transit-like events are naturally weighted much higher 
with this method. The fixed-transit duration allows a pri- 
mary sweep on all stars, which would be too computation- 
ally expensive if we varied the duration. A 1.5 Rj up planet 
with a one day period would create a 1.3-2.0 h transit du- 
ration, for stars of spectral type M5-K0. The same planet 
with a seven day period would create a 2.5-3.8 h duration. 
We have found that as long as the observed duration does 
not differ by a factor of two from the fixed duration, our 
algorithm can identify the transit. 

From this first period sweep, we compute the transit 
signal-to-noise for each star. The transit S/N is calculated 
from the fit of the data to a constant light curve as compared 
to a transit light curve. 

Following the first sweep, stars with a significantly bet- 
ter transit fit (~400 stars, S/N > 8.0) are subjected to an- 
other period and duration sweep, which refines the possi- 
ble transit parameters. Finally, the stars are then analyzed 
individually (in folded form and unfolded) to consider the 
possibility of a transit. Stars which have single faint points 
are rejected, as well as suspicious transits which occur only 
on nights with known problems. 



4 DETECTION SIMULATIONS 

In order to estimate how many stars might yield planetary 
transit detections, we used Monte Carlo simulations on two 
CCDs to estimate how many transit-like events we could re- 
cover if every star had a hot- Jupiter sized planet. We ran the 
simulations on CCDs one and two, and found very similar 
results. We assume that the other CCDs will show similar 
results, because all CCDs have similar magnitude distribu- 
tions. 

We began by randomly assigning each star a planetary 
inclination, planetary period, and planetary transit epoch. 
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The inclinations were uniform in cosi, the random period 
was uniform in logp from 3 - 5.2 d, and the epoch of mid- 
transit was a random date between zero and the period. The 
planet was assumed to be 1.5 Rj up and the stellar radius 
was computed using the colour information for each of the 
stars using equation 10. We then tested each of the systems 
to determine if the inclination allowed for a transit to be 
observed, and we compared the transit timing for each of 
the stars with our actual timings of our observations to see if 
the simulated transits would occur during our observations. 
Finally, we injected the transit into the data set using a 
simple box transit: if an observation was taken during the 
planet crossing the limb of the star, then the brightness of 
the star was decreased by half the full transit depth; if it 
was taken during the full transit, then the magnitude would 
be offset by the amount computed for a star that size being 
eclipsed by a 1.5 Rj up planet. 

After running this simulation on ~12,500 stars that 
were recorded on CCD one, we found ~720 stars (5.8%) that 
should have transits observable based on inclination and 
eclipse timing. Similarly, on CCD two we ran 14,000 stars 
and found ~800 (5.7%) that would transit. We then inserted 
these injected transits into our data set and loaded them into 
OPTPHOT, our transit search algorithm. We searched over 3 - 
5.2 d periods for 3 h transits. We were able to recover ~370 of 
the ~1520 stars with known transits (~25%). However, this 
does not suggest that our algorithm is missing well-defined 
transits. All stars were given a planet, and over ~55% of 
our stars are magnitude 21 or fainter, with an average preci- 
sion of 0.05 magnitudes. This precision at faint magnitudes 
prevents the detection of transits that would only produce 
shallow dips, especially since it would require many transits 
during our observing windows, an unlikely event. 

We are able to see a distinct differentiation between 
stars with injected transits and normal observed stars in 
our transit search. Figure 5 shows that the stars with an in- 
jected transit rise significantly above the stars without such 
a transit. This makes us confident that we would be able to 
find well defined transits in our brighter stars. 



5 RESULTS 

Presented in this section are the results from the observa- 
tions of NGC 6940. Similar results for NG C 7 789 or NGC 
6819 c an be found in lBramich et al] i2004l) and lStreet et alj 
(2003), respectively. 

5.1 Multiple Transit-like Events and Variable 
Stars 

Our transit search algorithm has discovered 14 stars in 
the field of NGC 6940 that have multiple short duration 
eclipses. Using the transit depth and stellar radii computed 
from their colour indices, we have determined a possible 
radius of each of the stellar companions. Every stellar 
companion is smaller than 35% the Sun's radius, and six 
are smaller than 25% of the Sun's radius. Folded lightcurves 
can be found in Fig. 6 while the parameters of each system 
are found in Table 2. The authors may be contacted for 
the complete data on each of the candidates to facilitate 
follow-up work. 



Star 6405 : Our eclipse depth of 8.9% appears to be too 
conservative, so the assigned companion radius of 2.1 Rj up 
is probably too small. However, the most damning feature 
of this eclipse (in terms of it being a planet) is the shallow 
secondary eclipse that occurs at half the orbital phase. This 
is definitely a binary system. 

Star 16016 : This is one of our best sampled transits, with 
4 transits observed. The noise amplitude of the lightcurve is 
consistent with other 20th magnitude stars in our sample. 
The eclipse bottom does not look particularly sharp, though 
sparse time sampling could have hidden that feature. We 
have computed a companion radius of 2.1 Rj up - If possible, 
this star should be measured using RV. 
Star 9939 : This has a fairly sharp eclipse, though it is 
only well sampled on egress, suggesting a grazing binary 
star. It is also fairly deep, with nearly a 25% drop in 
magnitude. However, our colours suggest the parent is an 
M2 star, with a radius of a little under half a solar radii, 
giving the companion a radius of ~2.4 Rj up - 
Star 13652 : This 18% eclipse is not as sharp as some of 
our other obvious binary stars, though the faint magnitude 
has introduced enough noise to make it difficult to ascertain. 
The parent star is one of our brighter candidates, a K4. The 
estimated companion radius is 3.1 Rj up , though that is a 
lower limit, as our eclipse may be deeper than our model 
suggests. Thus, it is probably a star. 

Star 1068 : If it were indeed a planetary transit, the 
companion radius would be around 2.2 Rj up , orbiting 
the K5 parent star. However, the sharp eclipse suggests 
a grazing binary, though sparse time-sampling and few 
observed eclipses may have contributed to that perception. 
Star 1254 : We cannot really classify if the eclipse is 
sharp or round bottomed, due to few eclipses and sparse 
time-sampling, though we think this could be a grazing 
binary. If the eclipse is caused by a planetary companion to 
the K5 star, R c would be around 3 Rj up - 
Star 2133 : Out-of-eclipse variation suggests that perhaps 
this is a binary star. However, it is one of our faintest 
eclipses, with a 3.9% dip found with three observed eclipses. 
This would indicate a planetary companion of 1.5 Rj up , 
which is well within the range for hot Jupiters. We suggest 
a follow up study of this star. It is also one of the brightest 
stars in our sample at 17.4 mags, which makes it a good 
candidate for further research. 

Star 11807 : This faint star seems to have a somewhat 
sharp eclipse, suggesting a grazing binary. The eclipse depth 
of 11.4% may be too conservative, so the computed value 
of the companion at 3 Rj up is probably too small. 
Star 13180 : This star exhibits some significant out-of- 
eclipse sinusoidal variation. The sinusoidal period (4.46 
days) appears to be slightly but significantly out of sync 
with the eclipse period (4.04 days). The variation may be 
star spot activity on the star. A sharp eclipse suggests that 
this is a binary star grazing its companion. 
Star 6716 :Sparse time-sampling prevents us from defini- 
tively saying this eclipse has a sharp bottom, but it appears 
so, suggesting a grazing binary star. The model suggestion 
of a 12.8% dip is conservative, so the computed companion 
size of 2.7 Rj up is a minimum. 

Star 7350 : This somewhat deep transit could be sharp 
bottomed, but the time-sampling is too sparse to say for 
sure. Since the parent is a relatively bright K6 star, we 
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Figure 5. Results of transit injected on CCD one. The left panel shows the results of the ~700 stars which had injected transits. The 
bulge indicates the easily determined transit signals. The right panel includes all stars on CCD one, including the original stars and stars 
with an injected transit. Very few stars achieve our cutoff of a S/N of 8. 



Table 2. System parameters of stars that show multiple transit-like events. Non- integer values of Nt mean that we observed partial 
eclipses. 
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26.3% 


3.3 


0.689(6) 


0.353(4) 


3.45 


57.568 


2 


20 h 34 m 54?87 


+28°04'13'.'61 



Table 3. System parameters of stars that show single transit-like events 

Star R mag R-I 5m 5t R* R c epoch RA Dec. 

(mag) (mag) (mag) (h) (Rq) (Rq) (HJD- 2451300) (J2000.0) (J2000.0) 

CCD 1 

1995 19.786(8) 0.956(21) 20.5% 2.8 0.581(7) 0.263(12) 51.653 20 h 33 m 40?32 +28°10'09'/67 
CCD 2 
CCD 3 

11284 19.190(7) 0.631(15) 10.6% 3.2 0.719(8) 0.234(19) 91.573 20 h 34 m 31?69 +28°30'i5('22 
CCD 4 

1533 19.706(8) 0.626(21) 27.4% 3.5 0.722(11) 0.378(16) 88.506 20 h 33 m 36!75 +27°59'44('47 

2510 19.693(10) 0.771(42) 15.7% 3.8 0.652(18) 0.258(29) 59.603 20 h 33 m 42!35 +28° 01' 29'.' 66 
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Figure 7. Unfolded lightcurves from single transit events. The vertical lines delineate the nights of observation. 



have very little scatter in our data points, and the model 
fits relatively well. 

Star 8837 : This sharp eclipse has some scatter out of the 
primary eclipse, and could have a secondary eclipse that we 
have not yet found. Also, the companion is computed to be 
larger than 3 Rj up , so is probably another star. 
Star 12930 : Though this is one of the faintest stars in 
our list of candidates, we can find the periodicity because 
we have luckily observed four transits. However, the scatter 
does prevent us from saying if the eclipse is sharp or round 
bottomed. 

Star 15028 : Our models fit this eclipse exceptionally well, 
but it is fairly deep at 25%, and suggests a companion 
radius of 3.6 Rj up - Further observations would be necessary 
to determine the shape of the eclipse. 



5.2 Single transit events 

We have also discovered several single low amplitude 
transit-like events. We are unable to estimate a period for 
these events, but we can use the colour indices to compute 
the radii of the stars and the companions. Complete 
lightcurves are in Fig. 7 and the parameters are found in 
Table 3. The authors may be contacted for the complete 
data on each of the candidates. 

Star 1995 : Our models fit this eclipse well within the 
limited time-sampling, but it is fairly deep. However, our 
colours indicate a late type star with a radius of 0.581 Rq 
which suggests a companion radius of 2.7 Rj up - 
Star 11284 : Again, sparse time sampling makes it difficult 
to characterize the shape of the eclipse. Our computations 
suggest a companion size of 2.4 Rj up - 

Star 1533 : This is a fairly faint star, at nearly 20th 
magnitude, so there is some amount of scatter in our data 



points. However, the late type star (K5) can produce this 
fairly deep eclipse with companion 3.8 Rj up , which is a bit 
large for a planet, and is probably a star. At other points 
in the data, there could be secondary eclipses that arc 
unresolved with our limited time-sampling, so this could be 
a faint binary. 

Star 2510 : This single transit eclipse could be sharp- 
bottomed, and scatter could obscure secondary eclipses. 
However, we have computed a companion radius of 2.6 Rj up - 



5.3 System Models: Checking Transit Duration 

Using the stellar parameters in Table 2, we attempted to 
compare our measured transit duration with a computed 
transit duration, based on the size of the star (derived from 
the colour index), the size of the companion (based on the 
measured transit depth), and the period. Table 4 reports 
these values of transit durations and the ratio between the 
two. This is not an absolutely rigorous check (because a 
missed transit could give us a spurious period), but it does 
give us some idea as to if the system which we describe is 
actually a possibility. 

We find that four of our stars, 6405, 1068, 1254 and 
13180, have observed values within ~20% of the computed 
value of the transit duration. We don't feel that this is an en- 
dorsement of these candidates as planets (indeed, we know 
6405 to be a binary system), but we feel it does probably 
eliminate the other systems from being planetary systems. 
Further, each of the four systems have companions com- 
puted to be between 2.0-3.0 Rj up , too large to be considered 
planets. 
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Table 4. Computed versus Observed Transit Duration 





n !(f\ 

P W 


At (h) 


At (h) 




6405 


1.42 


1.1 


1.4 


0.79 


16016 


2.17 


1.5 


3.0 


0.50 


9939 


2.20 


1.3 


2.1 


0.62 


13652 


2.22 


1.6 


4.0 


0.40 


1068 


7.14 


3.2 


3.8 


0.84 


1254 


4.90 


2.7 


2.4 


1.13 


2133 


3.74 


2.0 


3.4 


0.59 


11807 


5.82 


3.3 


2.6 


1.27 


13180 


4.04 


2.2 


2.1 


0.92 


6716 


3.65 


2.2 


1.6 


0.48 


7350 


1.77 


1.3 


3.0 


0.43 


8837 


3.54 


2.2 


3.5 


0.63 


12930 


2.67 


1.8 


3.2 


0.56 


15028 


3.45 


2.2 


3.3 


0.67 



5.4 Modeling the Planet Catch 

Our simulations suggest that if all our stars had a hot 
Jupiter, ~5.7% of stars would show an eclipse. That is, 
the planet's orbit and orbital inclination would allow us to 
record that transit with our observation regime. Our tran- 
sit searching algorithm has shown that it can find ~25% 
of these transits, if they are randomly distributed over the 
magnitude r anges we have in ou r data set. Finally, recent 
research by iFischer. et al] J2004) has quantified the rela- 
tionship between metallicity and planet frequency, allowing 
us to quantify how many planets we would expect in our 
sample, if it mimics the solar neighborhood. 

The Besangon model supplies us with metallicities for 
each star in the model, specific to our galactic coordinates. 
We use these metallicities because we are looking at field 
stars in the direction of NGC 6 940, instead of memb ers of 
that cluster. We have used the IFischer. et al] J2004) data 
to estimate the probability that each star in our model has 
a planet, based on its metallicity. We estimate that if the 
same planet abundance holds in the direction of NGC 6940 
as in the solar neighborhood, then ~ 2800 of our stars have 
Doppler-detectable planets around them, with periods up 
to three years. Our observation regime and the inclination 
of the system allow us to see the transit of 5.7% of those 
systems, or 160 stars. Further, if the periods of extrasolar 
planets are assumed to be uniform over log space, then be- 
cause we are only looking for hot Jupiters, and not planets 
with periods up to 3 years, we will only see 18% of the stars 
with planets, or ~29 transits. Finally, our transit detecting 
algorithm, when tested on all stars in our data set, was able 
to find 25% of transits, or seven transits. 

We have produced high precision photometry for 
~50000 stars in the direction of NGC 6940. If we use the 
Besangon model coupled with the IFischer. et al. (2004J) re- 
lationships, then we should find about seven hot Jupiters 
in our data set. However, we have found no convincing hot 
Jupiters. Nearly all of our 'candidates' are almost certainly 
grazing binary stars, though a few simply have too little 
data to define them. One of our stars (16016) has an eclipse 
that might be round bottomed, and the computed radius 
of the companion is 2.1 Rj up , which may be an M dwarf. 
Another (2133), which exhibits out of eclipse variation (sug- 
gesting a binary star) has a computed companion radius of 



1.5 Rjup, which is well within the range of hot Jupiter radii. 
We recommend these stars for further study. 

We can use Poisson statistics to estimate the signifi- 
cance of our null result, using: 

e~ a a x 

m = (ii) 

where a is our expected number of planets (seven) and x is 
our actual planet catch (zero). We use this to calculate that 
there is only a 9.12 x 10 -4 chance of finding zero planets 
when we expect seven. This gives us a 3.3er null significance. 

The lack of detections is surprising, even given the ex- 
pected metallicity distribution in the stellar field population 
we surveyed. The main systematic difference between the 
population stud ied here and the sol ar neighborhood sam- 
ples studied by IFischer. et al] i2004h is that our stars are 
predominantly late K or M dwarfs of 0.7 Rq or smaller. 
Radial velocity surveys have only discovered two M dwarfs 
harbouring planets, but that could be an observational bias 
against M dwarfs, which are often too faint for RV studies. 
Our results point to a lower incidence of hot Jupiters among 
late K and M dwarfs than among F or G dwarfs, regard- 
less of the metallicity. lEndl et all <|2003T) have embarked on a 
study specifically aimed at finding if the formation history of 
M dwarfs prevents planetary companions, though they have 
not finished their surveys. Our results thus suggest that hot 
Jupiters are less common around M dwarfs, and the lack of 
planets is not an observational bias. 



6 CONCLUSIONS 

We have obtained high precision light curves for ~50,000 
stars in the direction of the open cluster NGC 6940 using 
differential image analysis. We have used Monte Carlo sim- 
ulations to estimate how many transiting planets we should 
expect to find, assuming planetary frequency of the solar 
neighborhood. We determined the sizes of the stars using 
colour inf ormation from our observations and calibrated 
stars from Stetson (2000). Using a matched filter algorithm, 
we have identified several stars that exhibit behavior similar 
to that which is produced by an extrasolar planet. However, 
most of our candidates exhibit secondary or sharp-bottomed 
eclipses, suggesting that the stars in question are binary 
stars and not stellar systems with hot Jupiters. 

We have been unable to find the number of stars with 
transiting planets we estimated we would find. This could be 
because we are looking at mostly late type K and M stars, 
instead of earlier type F and G stars. We have identified sev- 
eral candidates with multiple transit-like events, and some 
with single events, though none are unambiguously caused 
by planetary companions. 
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